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Quintessence field is a widely-studied candidate of dark energy. There is "tracker solution" in 
quintessence models, in which evolution of the field (j) at present times is not sensitive to its initial 
conditions. When the energy density of dark energy is neglectable (Sl^ ^ 1), evolution of the tracker 
solution can be well analysed from "tracker equation". In this paper, we try to study evolution of 
the quintessence field from "full tracker equation", which is valid for all spans of Q.^. We get stable 
fixed points of and Q,^ (noted as and fi^) from the "full tracker equation", i.e., and Q.^ 
will always approach and Q.^ respectively. Since and are analytic functions of (j>, analytic 
relation of ~ Q.^ can be obtained, which is a good approximation for the ^ Sl^ relation and 
can be obtained for the most type of quintessence potentials. By using this approximation, we find 
that inequalities < and il^ < il^ are statisfied if the (or w^) is decreasing with time. In 
this way, the potential U{(j)) can be constrained directly from observations, by no need of solving 
the equations of motion numerically. 
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I. INTRODUCTION 

Present astronomical observations require the exis- 
tence of dark energy, a significant component of the uni- 
verse with a negative pressure 5|-[9| . Though it has been 
more than ten years since its discovery, one is yet to tell 
what the dark energy is. We are still analyzing proper- 
ties of dark energy from observational data and seeking 
suitable candidates. Most properties of dark energy de- 
pend on two parameters: the equation of state Wde and 
the fractional energy density fide- Once the Wde ^ ^de 
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relation is obtained, we know almost all we need. At 
present, it is still not possible to constrain the evolution 
of dark energy from observations jl0l - ll2| . There are only 



definite constraints of present values of Wde and fide from 
observations: w'"^^ is rather close to —1 and fl^^] is dom- 



inating (about 70%) 
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15| . More constraints on Wde 



and fide will be forthcoming from future observations, 
to get the evolution of the Wde ^ fide relation from the 
observations, more theoretical efforts should be made. 

At present, the most economical candidate of dark en- 
ergy is still the cosmological constant A, whose equa- 
tion of state w\ = —1. There is only a free parame- 
ter r^A in the flat ACDM model. But it suffers from 



several problems, such as the coincidence problem and 
the fine tuning problem. Another well studied candi- 
date is the quintessence 0, a slowly rolling scalar field, 
analogous to the inflaton. Its equation of state is = 
(02/2 - ?7)/(0V2 + U) so one has ~l < < 1. In 
quintessence models, the coincidence problem and the 



14| . For example, 



fine tuning problem can be alleviated 
there are tracker solutions for certain type of quintessence 
models, in which the evolution of (j) todayis not sensitive 
to its initial conditions at early times l|. The coinci- 
dence problem thus becomes less severe. But it is diffi- 
cult to find quintessence models with analytic solutions of 
equation-of-motion, due to the existence of background 
matters (dark matter, baryon and radiations). To study 
evolutions of quintessence models and to be compared 
with observations, one usually has to solve the equations 
numerically. There are efforts to find analytic approx- 
imations for solutions of equations of motions, such as 
16l | which gives a first order approximation solution for 
inverse power law potentials. 

In this paper, we will try to approximate the ~ fi^ 
relation at the recent fl^ dominating period in a semi- 
analytic way. To make sure that the evolution of (f) 
at present only depends on U{(j)), we assume there was 
tracking solution at early times. In [l|], conditions for the 
existence of tracker solution was given by the "tracker 
equation" , which is a differential equation for . But 
this "tracker equation" are only valid as $7^ <C 1. For 
our purpose, we need a full tracker equation that is valid 
for all Cl(h without conditions attached. Such an equation 



has been obtained 2tJ1 and will be used here to study 
evolutions of quintessence models. 

The paper is organized as follows. In section II, 
we introduce two new functions and Q,/, which are 
fixed points of the full tracker equation. Assuming that 
r = U"U/U'^ and e = {U' /UY/2 are nearly constant, 
we find that the fixed points are stable for and fi^ 
if r > 1. If r and e do not evolve extremely fast, 
the relation of ^ fJ^ will always approach to that 
of ^ Q^. In section III we show comparisons be- 
tween {{^^{(j)), ri0(0)} and {it;0((/)), fl^{(f>)} numerically 
for several typical quintessence models. The relation of 
~ is shown to be a good approximation for the 
^ flfj, relation. In section IV we show how to con- 
strain U{(t>) directly from observational conditions on 
and Q,f, through Wij, and fi^. Observational conditions are 
converted to simple inequalities for U{(t>). We conclude 
in section V with discussions. 



II. GET THE APPROXIMATION OF ~ 
RELATIONS 



The equations of motion for quintessence field are 



'<j) + + [/' = 

H2^(f)2 = l(p,„+p, + i<^2^f/) 



(1) 



from which one gets the equations for and 



3(1_Kj0) i 2 



r- 1 



(2) 



2(1 +UJ4,) (1 + cj^) (6 + .t)2 2{l + uj^) 6 + i 



{l+uj^){Q + x) 



(3) 



where 



1 



din a; 



1 — W0 2U' din a' dlna^ 

and a is the expansion factor. We have assumed a flat 



universe (fifc + fl^ — 1) and set Mpi = 1/\/8ttG — 1. 
The subscript b represents the dominating background 
matter. As fi^ ^ 1 at early times, Eq.® reduces to the 
"tracker equation" in l| . At the recent acceleration era, 
51^ is dominating and can not be neglected. One must 
use the full tracker equation Eq.Q. Note also Wb = 
in this case. In this paper, we assume that there was a 
long enough tracking period at early times, so that the 
evolution of the field at present depends only on U{4)). 
Eliminating il^ in Eq.Q by using Eq.(l2]), one gets: 
uj^ 2 X 



r- 1 



2(1 (1+w^) (6 + i;)2 
1 — X ^(p /a • \ 



(4) 



2(l + w^)6 + i 8e' ' 
For constant e and F, the fixed point (also called critical 
point) of Eq.Q (obtained by setting x ~ and i = 0): 



UJ4, 



^ (-3 - 2e + 4eF - ^/{3-2e + Ae^y 



24e 



is stable only if 



F > 



5 + 3^4 



(5) 



(6) 



6 + 2^4, 

where the $7^ value of the fix point is obtained from 
Eq.dS]) and ^ (also setting i = 0): 



n^^l-(s-2e + 4eF - a/(3 - 2e + 4eF)2 - 24e) (7) 



When Eq.® is satisfied, fl^, is also stable. In this 
case, aj0 and fl^, will always approach Q^, and fi^ re- 
spectively. In this paper, we will only study the case of 



F > 1 {i.e., W4, < wt), so Eq.® is guaranteed for all 
spans of il^. 

F and e generally are not constants, as they are func- 
tions of U{(t)). The above results are still valid if the 
evolution of w^, is not extremely fast, which can be sat- 
isfied in the most quintessence models. In this case, Wcf> 
and will keep on chasing the dynamic and O^. 

Giving the form of U{(l)) of a quintessence model, one 
gets parametric functions Wcj,{(t>) and ^cj,{4>) from Eq.Q 
and ([7]), and thus the analytic relation of Wfj, ^ fi^. For 
certain models, there are simple and explicit relations of 
W4, ^ flfj,. For example, for power law potentials U = 
{n > 0) one has: 

1 



W0 



(8) 



1 + ?i(l - l]0)/2 
The w^i ^ relation is a good approximation for that 
of ~ ri^, as the evolution of Wij, ^ Q^, will approach 
that of Wfj, ^ flfj,. We will show this in the next section. 
In this way, evolutions of quintessence models can be 
studied directly from U{(j)). 

III. COMPARED WITH NUMERICAL RESULTS 

In this section we will show that and fi^ are good 
approximations for w^, and fJ^, and so is the w^, ~ fl^ 
relation for that of ^ fi^. We have checked it for a 
variety type of quintessence potentials, and typical exam- 
ples are shown in Fig. [T]and Fig. [2j The accuracy of this 
approximation is precise enough to study the evolution 
properties of quintessence models, especially the models 
that are favored by present observations. As must 



decrease from its tracking value (close to Wb) to present 
value (close to —1), we will only study models in which 
decreases monotonously (i < 0). 
At first we estimate differences between and and 
between 17^ and fJ^. The Eq.(l2]) can be rewritten as: 



(i^)^(§^)-'=(i+i)^ 



(9) 



For a variety of quintessence models, we have seen nu- 
merically that (1 + w^)/{l + w^) and have the 
similar evolving forms as that of (1 + and 1 > 
> (1 + x/&f > (1 + w^)/{l + w^). Typical ex- 
amples are shown in Fig. [TJ If the evolution of (and 
w^) is slower, the value of x will be closer to 0, and the 
differences between w^, $7^ and their fixed points will be 
smaller. 

There is a lower bound x > Qw^/{1 — 2w^) given in 
2, l3|. As r^^/il^ is much closer to 1 compared with (1 -|- 
w^)/(l + w^), one gets a upper bound for the deviation 
A of from by setting Vl^/Vt^ ~ 1 in Eq.(|9l): 



[2 - iw^){l + w^) 



(10) 



' (1 - 2w^Y 

which is rather small when is close to —1, as shown 
in Fig. [TJ Present observations indicate that is rather 
close to —1 at low redshift. For most models A is much 
smaller than this bound as is not so close to —1, as 
shown in Fig. [TJ The deviation of fi^ from fJ^ is also 
small. 

The ~ ri^ relation thus is a good approximation for 
the ~ ^Ifj, relation. Several examples are shown in Fig. 



[2l At the early tracking era, Vl^ «1 and the relation of 
~ fi^ is almost the same as that of ^ D,^. When 
il^ becomes unnegligible, the curve of tw^ ^ i}^ will begin 
to get away from that of fl^ in the w — space. 

The curve of ~ fl^ will chase after that of ^ il^. 
Normally will tend to —1 and will tend to 1 at 
last, and the two curves will be close to each other once 
again. 

Empirically, we have also found a better approximation 
for the relation of D,^ on the basis of and fi^: 

_ (2 + - - 1 ~ - 

w^^w^^ [l + w^), n^^ii^ 

5 — Sflcj, 

(11) 

The curve of uJ^ ~ fi^ is much closer to that of ti;^ ^ Q^, 
as shown in Fig. [2J 



IV. CONSTRAIN QUINTESSENCE 
POTENTIALS 



In the above, we have obtained approximations 
and Q,f, for and which are analytic functions of 
U{4>). We will show how to constrain U{4>) directly from 
observational results on Wde and Qde through and 
Present data seems to indicate that < —0.8 

and 0.7 < < 0.8 
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15|. As more conditions on 



dark energy to be obtained in future observations, more 
quintessence models can be checked with directly by us- 
ing our method. 

At the early tracking era was close to — 1/3 



17 



18| , and present w^^^ is very close to — 1 . Taking this 
for guidance, here we consider only quintessence models 
in which (and {D^) keeps on decreasing monotonously 



(i < 0). This is guaranteed if U{(j)) satisfies the equation: 

(12) 



dln(r-l) 3, 1 , 

^ < — (1 

dlnU 2e^ 2r-l' 



In this case, one finds the following inequalities 



(13) 



if the evolution of is not extremely fast. Intuitively, 
the curve of ~ fi^ is always on the up side of that of 
~ 51^ in the w — Q space, as shown in Fig. [5) 
With the help of inequalities ([T^. J7((/)) can be con- 
strained directly from conditions on {wde^^de)- Take 



(14) 



for illustration [19]. Since vu^ decreases monotonously as 
ilcj, increases, we thus have 'w^{Q.^ = 0.8) < —0.8. This 
inequality can be converted to: 



r(e = 3/8) > 7/5 



which is a necessary condition for inequalities (1141) . 



(15) 



If is too close to — 1 , it will be difficult to distinguish 



3- 



quintessence models from the cosmological constant 
Take > —0.95 for illustration. It is then easy to see 
that Wci){^(t, — 0.7) > —0.95 is a sufficient condition for 



{wf^ > -0.95, > 0.7). Equivalently, 



r(e = 3/28) < 77/20 



(16) 



Listed in Table I are the constraints on parameters of 
typical potentials by Eq.([T5]) and p^ . 

We note that for certain potentials H.^ will tend to a 
maximum Qmax smaller than 1 at last, such as U{4>) — 
C/oe'^'/V0" (n > 0,0 > 0) [20]. These potentials always 



have a positive minimum Umin at a finite cf). According 
to Eq. ([7]) , as the potential rolls to Umi7i , ?7 = er will tend 
to a nonzero minimum rjmin with F — > oo and e — )• 0. In 
this case, Eq. (fT5|) and (fT6|) are still valid though il^max 
may be smaller than 0.7. 

V. DISCUSSIONS 

We have gotten stable fixed points and from 
the full tracker equation, and shown that they are good 
approximations for and Qrf, even in the fJ^ dominat- 
ing period, and are analytic hmctions of U{(f)). 
The relation of ^ fi^ thus is gotten from the para- 
metric functions (^^{(f) and ((/)), which is also a good 
approximation to the relation of Wcf, ^ fJ^. 

Formally, functions of Wfj, and with respect to ex- 
pansion factor a can also be obtained. Substituting 
Eq. ([5]) , (I?!) into the equation 
dfl(h 



(17) 



din a 

one gets the function of the field with respect to a upon 
integration. For example, for U = Un/(f>'^ (</> > 0) one 
has: 

Vu. 



3V5 



(18) 



where we have set present ri^"-* — 0.7 and ao — 1. Sub- 
stituting 0(a) into Eq.® and Eq.© one gets: 

= — (Vl20a3 + 49a6-7a3) 



1 



7 



(19) 



2 2Vl20a-3 + 49 
For most potentials, it is not easy to get explicit functions 
of 0(a), Wcf,{a) and flci>{a). 
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TABLE I: Constraints of typical potentials of quintessence 



U{(j>) (n > 0,0 > 0) 



— l(LL\^ 

— 'A U I 



V = 



&l ^ 5 2SJ ^ 20 



2p 



1 + -^ 

1 + ^ 

' n 

1 4- 

J- + (n-0^)2 



n > 



n > 1 



The critical points and can also be used to con- 
strain the potential of quintessence directly from obser- 
vational conditions on [wdeT^de)- We have adopted two 
conditions on present [w''^^ , ^'"^^ ) for illustration. Fur- 
ther astronomical observations will yield more properties 
of dark energy. It may give conditions on {wdeT^de) at 
other redshifts, or even the exact shape of the Wde ~ ^de 
relation. In that case, our method can be still usable 
to constrain the potential and study the properties of 
the quintessence models that are fit with observations 
directly. 

In this paper, we have only studied the case that 
(and w^) keeps on decreasing monotonously, from 
which the inequality ()13p is obtained. In fact, there are 
quintessence models in which is increasing at present. 
One example is the case with [/(</)) = C/o(e~'^/^ + e^^""^) 
2l| . In this type of models, will decrease to a min- 
imum close to —1 and then begin to increase. So the 
boundary for thawing and freezing fields in [19] will be 
crossed, as shown in Fig. |3l It can be shown that when 



dln(r-l) 3 
dhiU Ye 



W(j, will be increasing, so will be w^. It requires a rapid 
decrease of F. As F at early times must be close to 1 
to get enough tracking, usually there is a rapid increase 
of F at recent times. In this case the lower bound w' > 
— (1 — w){\ + w) for quintessence models 0] may be 
crossed too. It is because w — {wb — 2F -I- 2)/(2F — 1) 
will no longer be larger than if the increase of F is 
too fast. It can be seen in Fig. |3]that the line oi w' ^ w 
with the double exponential potential is ver y cl ose to the 



strict lower bound w' > 3w(l + w) given in 19|. For this 
type of potential, as and are increasing, there is 
an inequality similar to 1)13^ : 

W0 > ri^ > (20) 

This inequality can be used to constrain U{4i) from con- 
ditions on (wdc^de)- The methods used in this paper 
can also be extended to Phantom and K-essence models. 
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FIG. 1; Evolution of (1 + i/6)^ (solid lines), (1 + W4,)/(l + kj^) (dashed lines) and il^/fi^ (dotted lines) with respect to Q,^. 
The potentials: I. U = Uoe^^'^; II. U = Uo/^P; III. U = Uo/cfP'^. The last figure shows the deviation A of from for these 
models. 




FIG. 2: The ~ fl^ relation (solid lines), the ~ Cl^ relation (dashed lines) and the ~ relation (dotted linos) in the 
w — D, space. The potentials: I. U = f/oe^^"*; II. U = Uo/ff)^; III. U = (7o// ^ The last fi eure shows the differences between 
the relation of w,), ~ f?^ and that of Hj^, ~ iltt, for these models ( 6 = — [w(0) — w(f2 = iV)]/w{Q) ). 
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dw/dlna 




-0.95 -0.9 -0.85 -0.8 



FIG. 3: The curve of the quintessence model with U{(j}) = 
Uo{e~'^^^ + e"^"*^) in the w' — w phase space. This curve 
crosses the the boundary for thawing and freezing fields [l^ 
and the lower bound w' = —(1 — w){l + w) in 



